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We make a simple ansatz for the supersymmetric lepton-number violating Yukawa couplings, by 
relating them to the corresponding Higgs Yukawa couplings. This reduces the free B3 parameters 
from 36 to 6. We fit these parameters to solve the solar and atmospheric neutrino anomalies in 
terms of neutrino oscillations. The resulting couplings are consistent with the stringent low-energy 
bounds. We investigate the resulting LHC collider signals for a stau LSP scenario. 



I. INTRODUCTION 

The first experimental evidence for physics beyond the 
Standard Model (SM) has been found in the neutrino 
sector [l| . The solar and atmospheric neutrino anoma- 
lies are best explained in terms of oscillating mas- 
sive neutrinos |2|], as opposed to for example lepton 
flavour- violating interactions Assuming massive 
neutrinos, based on a three neutrino fit including the 
recent MINOS 3 and the SK-II atmospheric data [H], 
the corresponding neutrino mass and mixing parame- 
ters at lcr (3a) C.L. are [E 0] 



Ami 



\Am 2 31 \ 



7 o +0.27 
'• y -0.28 



+1.1 
-0.89 



x l(T 5 eV 



= 2.6 + 0.2(0.6) x 10- 6 eV z , 



33.7±1.3(1|;1 



43.3 



+4.3 
-3.8 



-9.8 
-8.8 







+5.2 / +11.5 



-0.0 



-0.0 



(1) 
(2) 
(3) 
(4) 
(5) 



The angles are given in degrees. The most widely dis- 
cussed extensions of the SM to include massive neu- 
trinos involve the see-saw mechanism Q- These re- 
quire right-handed neutrinos, as well as a new, typ- 
ically very large Majorana mass-scale. The see-saw 
mechanism can also be incorporated into the minimal 
supersymmetric SM (MSSM) [|], now requiring right- 
handed neutrino superfields, as well as the additional 
high mass scale. 

However, within supersymmetry there is another, 
in our opinion simpler, possibility to include massive 
neutrinos, namely via renormalizable lepton-number 
violating terms, Wy i7 in the superpotential (Toj . 
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h£ QiH d Dj 
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where we have employed the conventional notation for 
the superfields [Tl(. For later use, we have also in- 
cluded the superpotential terms involving the down- 
like Higgs superfield. The terms in Wy Ji violate R- 
parity (a Z 2 -symmetry) as well as proton hexality 
[121 [l3l | (a Zg -symmetry), but conserve baryon tri- 
ality (£?3, a Zs-symmetry, sometimes also mislead- 
ingly called baryon parity) [IE EE EH- The Majorana 
neutrino masses are generated via tree-level mixing 
with the neutralinos, as well as via radiative correc- 
tions [IE EE EE El, EE EH, H2- There is an implicit 
see-saw mechanism in the ncutralino-neutrino sector: 
nf/M 1 /2, but with a much smaller hierarchy of mass 
scales. Furthermore, no new fields or mass scales are 
required. 

Within a baryon triality supergravity model the 
largest neutrino mass is naturally small [ll| . For uni- 
versal soft breaking terms, the mixing, «j, with the 
neutralinos is zero at the unification scale. It is subse- 
quently generated at the order of a few MeV via renor- 
malization group equations. It is thus proportional to 
the product of a (small) down-like Higgs Yukawa cou- 
pling (for example of the bottom quark or the tau lep- 
ton), a (small) baryon triality coupling and the Higgs 
mixing parameter /i [H EE EE EHl ■ The lighter neu- 
trino masses are generated via radiative corrections, 
and are naturally further suppressed. 

There are 9 + 27 + 3 = 39 lepton-number violat- 
ing (complex) parameters in the superpotential Wp- 
There are also 39 corresponding soft-supersymmctry 
breaking parameters, which in principle arc indepen- 
dent, but are usually related to those of Wp via univer- 
sal soft-supersymmetry breaking (26[. In a top-down 
approach, e.g. based on the Froggatt-Nielsen mech- 
anism [23], one can attempt to predict the order of 
magnitude of all superpotential parameters, i.e. W^ p 
together with the Higgs Yukawa couplings, based on 
a spontaneously broken gauge symmetry, Ref . 1281. 12911 
and references therein. See also Refs. [2l|, EE l3fj. l3l|. 

In this letter, we instead propose a baryon tri- 



ality model of neutrino masses, based on a simple 
phenomenological ansatz, which relates the Higgs su- 
perpotential parameters to those that violate lcpton- 
number. The justification for this is that the down- like 
Higgs doublet supcrficld and the lepton-doublet super- 
fields have identical Standard Model gauge quantum 
numbers. We make no assumption about the possible 
underlying theory at the unification scale. This ansatz 
dramatically reduces the number of free parameters. 
If experimentally confirmed it would give a clear in- 
dication on how to construct the more fundamental 
unified theory. 

In the literature there are other simple ansatze (32l . 
[33l . l34j , the most common and also the most similar 
to ours is pure bi-lincar lcpton- number violation, i.e. 
\jk = Kjk = 0, and Ki ^ 0. For this there is an 
extensive literature, see for example [13, EH, HH, HO, 
l35l . [HI and references therein. We discuss how our 
ansatz differs from the bi-linear case in Sect. Hill 

A special feature of baryon triality models for the 
neutrino masses, is that they lead to other observable 
effects at colliders and can thus be tested I34l [371 l38l 

H m El, E2, 0, H S HE Ezl H Si, M Mm- 

In the case of pure tri- linear couplings (k^ = 0), a fit 
to the neutrino data, Eqs. (HJ)-©, leads to values in 
the range X ijk , X' ijk ~ 10" 5 - 10" 4 [H, 0, M, 
These couplings are very small, in particular, too 
small for the resonant production of supcrsymmctric 
particles [57j • However, they do lead to the decay 
of the lightest supersymmetric particle in the detec- 
tor, possibly with a detached vertex. This model can 
be confirmed by measuring the branching ratios of the 
various lightest supersymmetric particle (LSP) decays 
and thereby measuring the couplings. However, sev- 
eral points have in my opinion been missed in the 
literature. In the case of a pure fit, i.e. not a model, 
it is possible to have larger couplings, which do not 
contribute to the neutrino masses, or which arc not 
required for the fit. In this case, the LSP decay which 
dominates the collider signals will be completely inde- 
pendent of the neutrino sector. Thus pure fit models 
can only be tested if the neutrino mass parameters 
dominate the B3 sector. We consider here a complete 
model, where the fit to the neutrino data fixes all the 
B 3 parameters. Second, it has hitherto been assumed, 
that the LSP is the lightest neutralino. We go b eyond 
this and also consider a scalar tau LSP (Til |58| . 

Our analysis is structured as follows: In Sect. [HI we 
present our model in detail. We then briefly review the 
neutrino masses in baryon triality models, Sect. IIVI 
In Sect. [V] we estimate the values of the free parame- 
ters which result in acceptable neutrino masses. In 
Sect. IVI1 we numerically evaluate the new param- 
eters in our model, such that the neutrino masses 
and mixing angles fall in the required experimental 
ranges, cf Eqs. HJ-©. In order to obtain at least 
two non-vanishing neutrino masses, we must violate 



at least two lepton numbers. This typically leads to 
significantl y s tricter bounds on the products of cou- 
plings [59l. l60l l6ll . |62| . In Sect. IVIIi we investigate, 
whether our model is consistent with these bounds. 
In Sect. IVIII1 we discuss possible future tests of the 
ansatz at colliders, in particular the LHC. In Sect. IIXI 
we conclude. 



II. SIMPLE B3-MODEL 

In the MSSM, the lcpton doublet supcrficlds Li and 
the down-type Higgs superfield Hd have the same 
gauge quantum numbers. They are distinguished 
through a discrete symmetry: lepton number. How- 
ever, in the case of baryon triality, lepton number is 
violated and not well defined. In the most general 
baryon triality superpotential with the MSSM super- 
fields, Hd and Li have exactly corresponding terms in 
the superpotential, as can be seen in Eqs. ([6]) and (J7J. 
We take this correspondence to motivate the following 
simple ansatz for the Yukawa coupling constants 
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Here, £i, are c- numbers. Eq. §8§ has the required 
form to maintain the anti-symmetry of the Xijk in 
the first two indices. The ansatz for the dimensionful 
mixing terms, Eq. (JTU)), is no simplification and we 
retain the free parameters. Given the ansatz of 

Eqs. ([5]), (0, and assuming we know the Higgs- Yukawa 
coupling constants (leading to the SM fermion mass 
matrices), then the 36 couplings Xijk, Kjk' are param- 
eterized in terms of the six numbers £j , t'j . 

Since Li and Hd have the same gauge quantum 
numbers, our ansatz in Eqs. ([5]), ^ is given in 
the SU(2) x U(l) citrreni-eigenstate basis. Thus 
when computing neutrino masses and comparing the 
required Yukawa coupling constants to low-energy 
bounds, we must rotate to the mass-eigenstate basis 
(63l . |64| . This requires a bi- unitary transformation in 
generation space. We shall denote the transformation 
of the left-handed and right-handed fermions (not su- 
perfields), respectively by 

e L = V e e' L , d L = V d d' L , u L = V u u' L , (11) 
eR = U e e' R , d R = U d d' R , u R = U u u' R , (12) 

where V e ,d,u, U e ,d,u are 3x3 matrices in generation 
space and the prime denotes the mass-cigenstates. We 
have combined the charged lepton and quark states 
into three-component vectors in generation space, e.g. 
bl = Hl,t~l). We rotate the sfermion part- 
ners by the same matrices in flavour space. By con- 
struction, these transformations diagonalizc the SM 
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Coupling 
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Numerical Value 
(tan /3 = 10) 
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-£ 3 \/2m ti /vd 


-6.1 • 1CT 3 •& 


A233 


t2\f2m T /vd 


1.0- 10 _1 -£ 2 



TABLE I: Predictions for the LLE in our ansatz as a func- 
tion of the free parameters £i . 



Yukawa coupling matrices 
■ (h E f ■ V e 



— diag (m e ,m M ,m r ) , (13) 

Vd 



Uj-(h D ) T -V d 



Uj ■ (h u f ■ v u 



V2 

— diag(m d ,ra s ,m b ) , (14) 

Vd 

— diag(m U) m c , m t ) , (15) 
v,. 



where the normalization of the Higgs vacuum expec- 
tation value is v = \/\v u \ 2 + \vd\ 2 = 246 GeV and 
tan/3 = v u /vd- 

In the following, we assume that the charged lepton 
mass- and weak-eigenstates are the same. The cor- 
responding charged lepton rotation matrices are then 
given by Ve. = U e = 11. Thus in our ansatz, in the 
leptonic sector the mixing takes place entirely in the 
neutrino sector, cf. the discussion in Refs. [6l|, [64| . 
Using Eqs. ([H]), and (fl"3"j) . the LLE couplings can then 
be expressed in terms of the lepton masses and the 
three parameters £i 



v d 



V2r, 



Vd 



Oik 



(16) 



where m e j = (m e ,m^,m T )\ Explicitly the couplings 
are given in Table [H as a function of the free parame- 
ters li. As an example, we have also given numerical 
coefficients in the case where tan/3 =10, which fixes 
Vd = 24.5 GeV. Overall, of course, all couplings are 
proportional to the free parameters li. 

We thus have some very specific predictions for the 
LLE couplings in our model. 
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Besides the vanishing couplings, we would thus expect 
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the couplings to satisfy the strict constraints 
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100 GeV 



A 122 < 3.6 -lO" 4 
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100 GeV 
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(19) 
(20) 
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100 GeV ' 

where we have implemented the low-energy bounds in 
[rJlj for A133, A233, A232 and inserted the PDG lepton 
masses 



65] 



Throughout we consider our model only at the weak 
scale. In principle it should be embedded in a unified 
model at the grand unified scale or above [28, 29]. In 
that case, the predictions in Eqs. (fT7|) . (fT8|) would be 
modified by rcnormalization group effects. In partic- 
ular the couplings in Eq. (|17[) would get non-zero con- 
tributions [251 ] . which however are extremely small, as 
they are proportional to the product of three non-zero 
LLE couplings. 

When expanding the LLE term in the Lagrangian 
into its mass-eigenstate components, we obtain (sum- 
mation over generation indices implied) 



LLE 



-A 



ijk e R 



k *u ic P L e j - \ ijk e{e k P L S 



-\ ijk v l L e k P L e> 



+ h.c. . (22) 
Next we consider the LQD term in the Lagrangian. 



Expanding out the SU(2) l doublet superfields, we ob- 
tain 



■&LQD — K 3 k N LiD L jD Rk - \' tjk E Ll U L ]D Rk 



(23) 



Rotating the quark superfields in the first term into 
the superfield basis where the quarks are in the mass- 
eigenstate and using Eqs. © and (fl"3"|) we obtain 

X' ijk N Li D Lj D Rk = K ]k [U\] rk [V d ] ]S N Ll D Ls 'D' Rr 
= l\ [U\{h D ) T V d ] rs N iL D' Ls D' Rr 

= ^^m d r6 rs N iL D' Ls D' Rr . (24) 

Vd 

For the second term in Eq. (f2"3")) . we obtain analogously 

Kjk E LiULjD Rk 

= H i [U d \h D ) T V d \ rt [V d W u ] ts E' Ll U Ls 'D' Rr 

V21: 
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md 



r [V% r El/Uls'D'j 



Rr 7 



(25) 



where Vckm = V u 'V d is the Cabibbo-Kobayashi- 
Maskawa matrix. Combining Eqs. (j24|) and (|25|) . we 
can then write the Lagrangian for the LQD interac- 
tions in the mass-eigenstate basis and expanded in 
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supcrfield components 1 

JSP, 



LQD 

~Kjk (j>ld j P L d k + d 3 L d k P L v l + £*v lc P L d 

+ KjkV r ) (e l L d k P L u r + u r L d k P L e* + d k R *e lc P L u' 
+ h.c. , (26) 

with the coupling defined by 
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We also introduce the notation 



Ojk 



A' 



ijk 



I 



,V2 



m d . 



[V] jk ■ 



(27) 



(28) 



Note that the (s)neutrino interactions are flavour di- 
agonal in the down-(s)quarks, whereas the charged 
(s)lepton interactions involve generation off-diagonal 
(s)quark interactions [64J. We have given an estimate 

of the couplings A' and A', modulo the li in Table HT1 
We have again assumed tan/3 = 10. Furthermore, 
we have taken the central PDG values for the quark 
masses: m d = 6MeV, m a = 103 McV, m b = 4.2 GcV, 
and the central values of the global PDG fit for the 
CKM matrix entries (2 significant figures) (6f| 



0.97 0.23 0.0040 
Vckm = I 0.23 0.97 0.042 
0.0081 0.042 1.0 



(29) 



As can be seen from Table HI1 we thus have also sim- 
ple predictions for the A' couplings in terms of quark 
masses and Vckm entries and independent of tan (3, 



ijk _ m d kVjk 
m d iVji 



A' 



A! 



0' 



a- 



Vjk 
Vik 



We can use Eqs. ([15]). ([27]) l[2"5jl to translate between 
our parameters £i, l\ and the B3 couplings, where 
in the latter case, care must be taken to include the 
CKM-mixing for the charged (s)lepton interactions. 

It is the purpose of this letter to investigate whether 
with this reduced freedom in the B 3 sector, we can still 
obtain neutrino masses and mixings, which are, first 
of all, consistent with Eqs. (U)-© and second, where 
the resulting coupling constants are consistent with 
the existing low-energy bounds. In Sect. lVlITl we then 
study possible observable consequences of the absolute 
values of the couplings as well as of the relative values. 



The primes, denoting the mass-eigenstates, arc omitted in 
the following. 
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6.0 • 10~ 3 


5.8 • 10" 3 
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1.0 • 10" 2 


(131) 





2.8 ■ 10" 6 


(132) 





2.5 • 10~ 4 


(133) 


0.24 


0.24 



TABLE II: Predictions for the LQD in our ansatz as a 
function of the free parameters . In the right column we 
have assumed tan/3 = 10. 



III. OTHER ANSATZE 

In Rcf. [H, [H, HI the hierarchy in the SM Higgs 
Yukawa couplings was taken to motivate a similar 
hierarchy in the LQD (and separately in the LLE 
couplings). The authors restrict themselves to the 
couplings A'33 and A.;33. We extend this interesting- 
work in several respects. We include the most recent 
neutrino data in our fit; Eqs. |l])-([-3). Furthermore, 
we include the CKM mixing in our ansatz, we thus 
have a prediction for the full range of the couplings. 
This is particularly important for the observable con- 
sequences of the model, i.e. the LSP decays. We also 
do combined fits including all the couplings, i.e. the 
A and the A' couplings and also the «;,, our Models I 
and II below in Sect IVH 

The most widely considered simple ansatz are B3 
models, where A^fe = A^- fc = and k 2 ; 7^ 0, of- 
ten denoted bi-linear R-parity violation. This clearly 
has only three free parameters compared to the six 
or nine, in our models below. In order to compare 
the two ansatze in more detail, we combine the fields 
C a = (Cq, Ci) = (H d , Li), where a = 0, ...,3 and 
i = 1,2,3. The bi-linear i?-parity violating superpo- 
tential is then 

W = h% dCoEj + hE QiCoDj + n C a H u + Ki LiH u , 



We can now make a field redefinition 
C^C' = RC, 



(31) 



(32) 



such that the bi-linear lepton-number violating terms 
are eliminated from the superpotential. The explicit 
form for R is given in Ref. [TJ, [3, [28| . We then obtain 
the superpotential 



W 



[R]ia[R]o/3 CaLpEj 



I J 



[R]o a QiC a Dj + p, CqH u . (33) 
The transformed parameters (denoted by a tilde) arc 
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then given by 
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(36) 


A ijk 


= hf k [R} 0i 


[R]oi 

[R]oo 


hf k 


(37) 



\ijk 



'Ik 



[R]ij[R]oi — [R]n[R\oj 



\R 



:]«(^[«]«)-[«]oi([fl]K^). (38) 



In the last equation we have included the parentheses 
to emphasize the sum over I. Note that not only are 
the tri- linear couplings Xijk, A^ fe generated, but also 
the Higgs Yukawa couplings are modified 



E,D 



■ E,D 



From Eq. (|3T|) we see that 



K jk = e'ihfk, with % 



m 

\R 



llll 



(39) 



(40) 



In the last equation, we have employed the explicit 
form of the matrix R. Next we would like to show 
that 



Xijk = Uhfk-ijhZ. (41) 
For this we insert h E from Eq. ([56]) and factor £ij 

X ijk = £ 3 [R]oo(hf k [R}u) - e t [R]oo(hf k [R}i,) 

+ti [R) 0i (hf k [R) l0 ) - £ % [R] 0j (hf k [R} 10 ) (42) 

If we now set 

\R\0i , e [R]oj 



\R] 



and 



on 



\R 



llll 



(43) 



the last two terms in Eq. (|4"2"|) cancel and the first two 
terms agree with Eq. (|38| . In particular, we see that 
B3 models with only bi-lincar terms arc a special case 
of our ansatz with L = £'-. 




FIG. 1: The radiative slepton-lepton and squark-quark contri- 
bution to the neutrino mass. 



in one massive neutrino at tree-level. The resulting 
mass matrix is perturbatively given by 2 [H, H3, H(| 



i%M^ cos 2 f3 



ft (m 2 z M^ sin 2/3 - M 1 M 2 fi) 



\l ! 



(44) 



Here mz and M\p denote the mass of the Z a 
gauge boson and the soft supersymmetry breaking, 
electroweak gaugino mass parameters, respectively. 
The photino mass is given by My = Mi cos 2 0w + 
Mi sin 2 8\y , where 9w is the electroweak mixing an- 

gie. \n\ = #+rr 



|Ki| 2 . We shall see below that 
0( 1 MeV 2 ) and thus to a high precision 
H [ll|, 0, [23[. For later convenience, we have 
introduced the constant C3 to summarize the depen- 
dence on the supersymmetric parameters. Within 
B3-supersymmctry, the other two neutrinos obtain 
masses through radiative corrections from both the 
bi-linear [l7l.Jl8l and the tri-linear terms in the super- 
potential [3Hfl|- I n the following, we focus on the 
radiative corrections due to the tri-linear terms, since 
we expect these to dominate for small Ki [TlT | , however 
realistic neutrino mass models based only on bi-linear 
terms have been constructed [13, EH, HH • 



There are two distinct radiative contributions to the 
neutrino masses from the tri-linear couplings for which 
the Feynman diagrams are shown in Fig. |T]). One is 
proportional to X' ikt X'j tk , where a squark and a quark 
propagate in the loop. A second is given by a slepton- 
lepton loop, and is proportional to Xi k eXj£ k . For the 
squark-quark loop, the bottom-sbottom contribution 
dominates [k = £ = 3), since for the down- like quark 
masses we have m 2 3> to 2 , to 2 ,. This results in the 
mass matrix 3 



{mi)i 



2NcX' i33 X' j33 
16tt 2 



mlAb 



M? 



(45) 



IV. B 3 NEUTRINO MASSES 

The general, lepton number violating superpotential 
in a i?3-modcl is given in Eq.([6|). As stated in the 
introduction, due to the Ki, the neutrinos mix with 
the Higgsino components of the neutralinos, resulting 



2 As in O, we assume here a rotation into the basis where 
the sneutrino vacuum expectation values vanish. In principle 
this requires a detailed minimization of the scalar potential 
as discussed for example in Ref. llll 1 1 71 . |21H . This is well 
beyond the scope of this paper. 

3 The extra factor of two arises from two distinct contributions 
to the mass which are equal for k = I [2di | . 
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where N c = 3 is the colour factor and 

x b hix b f M h 2 \ { ac\ 

f(x b ) = -, with x b = —f- . (46) 

x b -l \ M bJ 

M~ h , s — 1,2. denote the sbottom mass-eigcnstates, 
where Mg < Mg 2 . A b = A a b - /Stan/3, where A° b 
is the tri-linear soft-supersymmctry breaking bottom 
coupling. i,j = 1,2,3 are generation indices. In cal- 
culating explicit values for Eq. (|4"S"|). the mixing of the 
left and right handed bottom squarks is taken into ac- 
count but generation mixing is neglected due to the 
strict constraints from flavour changing neutral cur- 
rents. For later numerical estimates, we find that 

/(1) = 1, /(100)«4.65, (47) 

where the latter value corresponds to the fairly ex- 
treme value of Mr = 10 -Mr . Note, there are no neu- 
trino mass contributions proportional to m b m s . The 
subleading contribution is proportional to A^ 22 Aj 2 2 m s- 
This could in principle be dominant for A' i22 > 
(m b /m s )X' i33 ~ 45 ■ A^ 33 . However due to the rela- 
tion Eq. (|27|) . this is not possible in our ansatz, cf 
Table [HI 

The contribution from the slepton-lepton loop is 
analogously given by (N c = 1) 

K)„^^#. (48) 

Here m T is the tau mass and x T = (M^/M^) 2 , < 
Mf 2 are the stau masses. A T = A® — /2tan/3, where 
A® is the trilinear soft breaking term for the r. Since 
Xijk is antisymmetric in the first two indices, we must 
restrict the indices i, j = 1, 2 in Eq. (|4"5)) . For i = j = 
3 the leading contribution 4 is from the smuon-muon 
loop, proportional to A§ 22 m^. 

At any given energy scale, the mass parameters Ki 
in Eq. © can be rotated away (67[. Depending on 
the scale and mechanism of supersymmetry breaking, 
the corresponding soft supersymmetry breaking terms 
will then also vanish [ll|. When this occurs, we have 
pure tri-linear i^-models, which have been widely dis- 
cussed in the literature. The leading neutrino mass 
contributions must then arise solely from the above 
loop-diagrams. In the following, we shall thus discuss 
two models embedded in our ansatz. In Model I, we 
consider the case of pure tri-linear interactions, i.e. 



4 In principle, we could also have a contribution proportional 
to Ai22Aj22"i^ °c ^3 m fi/ v d' This can be large since £3 is a 
free parameter in our ansatz. However, in order to have a 
comparable contribution, we must have I3 Rj (m T /rnu^) 2 « 
300. Comparing with Tabic IJ we see that the resultin g A2 32 
typically violates the low-energy experimental bounds |6ll |. 



Ki = 0. The neutrino masses are then solely given 
through the combined loop contributions in Eqs. (|45l) 
and ([48]) . In Model II, we consider the more general 
case with 7^ 0, as well as A, A' 7^ 0. In this case, 
either the complete Majorana neutrino mass matrix is 
given by the sum of all three contributions, Eqs. (|44[>. 
([45]) and (j48]) or one of the loop contributions is ab- 
sent. In both models, the neutrino masses and mixing 
angles are then obtained upon diagonalization. 

In the following, we shall first estimate the resulting 
neutrino masses in both Models I and II and then de- 
termine the masses and mixing angles by numerically 
diagonalising the complete mass matrix. 

V. NEUTRINO MASSES FROM THE 
SIMPLE B3-MODEL 

We now insert our ansatz into the neutrino mass for- 
mula? above and give an estimate for the £i, £j. The 
tree- level contribution, Eq. (f44|) . is unchanged. The 
one-loop contributions, Eqs. (|4"5)) and (|48[) . are given 
by 

« ^^V^U^jCV, (50) 

where in the last equation i,j 7^ 3. For later conve- 
nience we have introduced the constants C , T . In or- 
der to estimate the order of magnitude of the ti, i'j, we 
shall assume hierarchical neutrino masses. In Model 
I, the neutrino masses are generated alone from the 
loop corrections and the mass matrix is then 

[M4 y . = tdjCr + U'^Cb , (51) 

which results in two massive neutrinos, since i, j 7^ 3 
for C T . We can obtain a third massive neutrino if we 
include the subleading term proportional to liljC^, 
for which there is a contribution for i,j = 3. 

In order to obtain an estimate, we set tan/3 = 10 
and assume A% = A% = M~ b - = m soft = 100 GeV, 
which results in f(x) ~ > 1. We then have C b ~ 130 
keV and C T = 4.1 keV. The heaviest neutrino is in 
agreement with the atmospheric neutrino anomaly, 
Eq. ©, for 

e w 3.5 -10" 3 , or £' «3.6-10~ 4 , (52) 
A l33 « 3.6 ■ 10~ 4 , or A^ 33 w 8.7 • 10" 5 . (53) 

Correspondingly, we can generate the mass required 
by the solar neutrino anomaly by the other term. We 
obtain m v w 0.01 eV, for 

I k, 1.6 -10 -3 , or t' wl.6-10" 4 , (54) 
A 433 « 1.6 -10" 4 , or A-33 w 3.9 • 10" 5 , (55) 
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For both the atmospheric and solar anomalies we have 
used Eqs. (fT6|) and (f27|) to translate back to the cor- 



responding values for A.^ and A' 



(3.3- 



In Model II, the largest neutrino mass is generated 
by the tree-level contribution. Taking the trace of 
Eq. (|44p. and assuming Mi = M 2 = fj, = m so f t , we 
have approximately 



"soft 



From Eq. ([T]), we then obtain [25 



= C(l MeV). 



(56) 



(57) 



With the same assumptions as in Model I, the ra- 
diative contributions then generate the solar neutrino 
mass for the values given in Eqs. (j54"|) , (|55[) . 



VI. 



NUMERICAL EVALUATION OF THE 
NEUTRINO MASSES 

A. General Outline 



Next we wish to determine more precisely the param- 
eters of our ansatz, i.e. the £i, i^K-i, by fitting them 
to the neutrino data. They in turn fix the B3 param- 
eters. In order to learn more about the importance of 
various parameters, we consider two cases. In Model 
I, we consider the case of pure tri-lincar terms. The 
free parameters arc 



Model I 



(58) 



In Model II, we include Kj 7^ 0; the respective free 
parameters are given by 



Model II 



(59) 



can also lead to an observable effect. If both the mass- 
less and the massive, second lightest neutrino have 
significant electron-neutrino admixtures, then the cor- 
responding Kuric plot will have a dip at the electron 
energy Q — 777,2, with 777,2 being the second lightest neu- 
trino mass. The maximal electron energy however will 
be Q, within experimental uncertainties [68) . Depend- 
ing on the parameter values this could be observable 
by the KATRIN experiment [11]. 

In Model II, we obtain three non-zero neutrino 
masses from the real symmetric mass matrix, 



M 



11 



( K\K 2 K1K3 

K 2 Kl K 2 K 2 
\ K 3 K 2 



Q, 



^ t'xt'i t'\t! 2 i'l&'s 
pi pi pi pi pi pi 

*-2*-l c 2^2 C 2 C 3 



^3^2 



+ 

\ 



(61) 



which depends on 8 independent parameters. 



In the numerical evaluation, the coefficients Cg, Cb 
and C T are determined by assuming a BC1 mass spec- 
trum [58|, which has a scalar tau LSP. The resulting 
neutrino mass matrix (M„)jj is diagonalized by the 
orthogonal rotation matrix V 



V T MLV = 




(62) 



where m\ < m 2 < 7773 and V is given by the standard 
parameterization 



V = 



(63) 



C12C13 
-S12C23 - C12S23S13 
S12S23 - C12C23S13 



S12C13 S13 
C12C23 - S12S23S13 S23C13 
-C12S23 - S12C23S13 C23C13 



Clearly, Model II is the most general case. 

The full neutrino mass matrix (Mj,)y is given by 
the sum of Eqs. (|44|) . (|49|) and (|50|) . where depending 
on the model, the coefficient of Cg can be zero. In the 
case of Model I, the real symmetric mass matrix has 
the form of 



M}, = 



(60) 



i 2 i x 



£ x i 2 
£' 2 £' 2 
^'3^2 



^1*3 



^3 



till hi* 
I £ 2 £i 




which is a function of five parameters. However, due 
to the simple form of the matrices, there are only two 
non- vanishing neutrino mass eigenvalues. This is suf- 
ficient to explain the atmospheric and solar neutrino 
anomalies. In principle, a masslcss lightest cigenstate 



with Cij = cos dij and s, j = sin 6ij . The complex Dirac 
phase (5i3 and the two Majorana phases ai/ 2 are omit- 
ted. The experimental ranges of the masses and the 
mixing parameters 6\ 2 and #23 are given in Eqs. {!])- 
In addition, the bound on the mixing angle #13 
from the CHOOZ experiment Eq. (JSJ) is taken into ac- 
count. We find it convenient in presenting our results 
instead of using the angles as in [1,0], to use tan 2 9i 2 , 
tan 2 023 and sin 2 #13. The corresponding ltr ranges 
are given by 

tan 2 (9 12 £ [0.403,0.490], (64) 
tan 2 6» 23 e [0.680,1.20], (65) 
sin 2 6» 13 < 0.0082. (66) 

In performing the fit, we randomly sample log 10 (^) 
and log 10 (^), thus guaranteeing that we explore the 
full hierarchy of couplings. We only consider couplings 
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&i, £'j > l ■ 10 7 , as smaller couplings have no effect 
on the neutrino observables. We furthermore require 

liij -C T < 0.1 cV (67) 
e'/j-Ct < 0.1 cV (68) 
KiKj ■ C s < 0.1 cV. (69) 

For given values of the parameters £i, £'p Kk, the re- 
sulting mass matrix is computed. The numerical diag- 
onalization of the mass matrix yields the mass eigen- 



values mi, m2, and the orthogonal transformation 
matrix V and thus the mixing angles (9 12 , #23 and #13. 
Afterwards, all experimental requirements, Eqs. JT])- 
arc applied on the mass eigenvalues and on the 
mixing angles. We delete models, which do not fall 
within the lcr ranges. Throughout we also assume a 
hierarchical mass spectrum. Our results are summa- 
rized in Tables IIIII and HV| which we discuss in the 
next section. 





£1 


k 


£'1 


£'2 


£'3 


tan 2 6x2 


tan 2 #23 


sin 2 <9i3 


Ami! 


Am| 3 


i\ max 


9.38 • 1(T 4 


1.69 ■ 10" 3 


-1.60- 10~ 7 


4.01 • 10 -4 


-5.22 ■ 10 -4 


0.48 


0.99 


0.0078 


0.081 


2.5 


i\ min 


7.10- 1(T 4 


-1.76 ■ 10~ 3 


9.27 ■ 10 -6 


4.16 ■ 10~ 4 


-5.09- 10~ 4 


0.42 


1.10 


0.0037 


0.077 


2.6 


£2 max 


-7.12 ■ ltr 4 


1.78 ■ 10~ 3 


-9.42 • 10~ 5 


-4.10 • 10~ 4 


5.04 • 10" 4 


0.42 


1.11 


0.0038 


0.081 


2.5 


£2 min 


9.23 ■ 1(T 4 


1.58 • 10" 3 


1.82 ■ 10~ 6 


3.83- 10~ 4 


-5.60- 10~ 4 


0.49 


0.72 


0.0056 


0.078 


2.7 


£1 max 


7.io- itr 4 


-1.76 • 10~ 3 


9.27 • nn s 


4.16 • 10 -4 


-5.09- 10" 4 


0.42 


1.10 


0.0037 


0.077 


2.6 


£1 min 


-9.14- itr 4 


-1.67- 10" 3 


1.00 ■ 10" 7 


4.11 • 10~ 4 


-5.35 • 10" 4 


0.47 


0.95 


0.0065 


0.077 


2.7 


£'2 max 


9.01 ■ 10~ 4 


-1.77 ■ 10~ 3 


1.10 ■ 10~ 7 


4.28 • 10 -4 


5.16- 10" 4 


0.43 


1.18 


0.0077 


0.078 


2.8 


£'2 min 


8.73 ■ 10" 4 


1.62 • 10 -3 


4.00 ■ 10~ 7 


3.61 • 10~ 4 


-5.53- 10~ 4 


0.41 


0.68 


0.0053 


0.080 


2.4 


£'3 max 


8.69 • 10~ 4 


1.61 ■ 10~ 3 


-5.50- 10~ 7 


3.82 • 10 -4 


5.71 • 10 -4 


0.42 


0.69 


0.0043 


0.078 


2.8 


£3 min 


-7.12 ■ 10~ 4 


1.78 ■ 10 -3 


-9.42 ■ 10~ 5 


-4.10 ■ 10~ 4 


5.04 • 10" 4 


0.42 


1.11 


0.0038 


0.081 


2.5 



TABLE III: Explicit solutions for Model I, where Ki = 0. We only show the values where one of the \£i\, \£'j \ takes on an 
extremal absolute value, i.e. the highest or lowest value obtained in our fits. Note that due to our sampling constraints, 
we did not probe values for the £i, £'j < 10~ 7 . In the five columns on the far right, we show the resulting neutrino 
parameters. The values for Am^ and Aml 3 are given units 10 -3 eV 2 . Comparing with Eqs. CD)-©, we can see which 
physical parameter is at its experimentally allowed limit. Thus for example in the first row, where £\ is maximal, tan 2 #12 
and Am 2 ,! are at the edge of their allowed values. Pushing £\ any higher would violate these constraints. 





£1 


£2 


£'1 


£'2 


£', 




K2 


K3 


£\ max 


1.08 • 10 -3 


-1.17- 10~ 3 


1.80 • 10 -6 


-4.39 • 10 -4 


4.57- 10" 4 


-0.128 


0.00551 


-3.40 


£1 min 


1.30- 10~ 7 


1.45 ■ 10~ 5 


-6.50- 10~ 7 


-4.30 ■ 10~ 4 


5.35 • 10~ 4 


1.93 


2.82 


1.09 


£2 max 


-8.15 ■ 10~ 4 


1.77 • 10 -3 


4.06 • 10 -5 


4.14- 10 -4 


-5.01 • 10" 4 


-0.159 


-0.541 


-0.103 


£2 min 


4.50- 10~ 7 


1.20 ■ 10~ 7 


1.42 • 10~ 4 


3.22 ■ 10~ 4 


4.10 ■ 10~ 7 


1.02 


-4.65 


6.11 


£\ max 


-8.80- 10~ 7 


1.34 • 10~ 3 


1.84 ■ 10~ 4 


2.28 ■ 10~ 4 


2.20 • 10 -6 


0.401 


-4.64 


6.01 


£[ min 


-4.65 ■ 10 -6 


9.62 • 10 -6 


2.20 • 10~ 7 


-5.28 ■ 10~ 4 


-4.04- 10~ 4 


-2.07 


-0.0320 


3.72 


£'2 max 


-1.50- 10~ 6 


-7.07- 10~ e 


-8.85 ■ 10~ 5 


5.31 ■ 10~ 4 


-4.10 ■ 10~ 4 


1.70 


0.100 


-3.81 


£'2 min 


-9.49 ■ 10 -4 


1.66 • 10 -3 


-9.70 • 10 -6 


1.30- 10~ 7 


1.10 • 10 -4 


0.219 


4.92 


5.98 


£3 max 


8.82 ■ 10~ 4 


1.64 ■ 10~ 3 


3.60- 10~ 7 


3.91 ■ 10~ 4 


5.55 • 10~ 4 


-0.457 


0.0604 


-0.218 


£'3 min 


4.50- 10~ 7 


1.20 ■ 10~ 7 


1.42 • 10 -4 


3.22 • 10 -4 


4.10 • 10~ 7 


1.02 


-4.65 


6.11 


ki max 


1.43 • 10~ 4 


-1.30 ■ 10~ 3 


-2.69 ■ 10~ 6 


4.47- 10~ 4 


-4.77- 10~ 4 


2.06 


-0.0125 


-3.03 


ki min 


-1.87- 10~ 6 


-2.30 • 10 -5 


-1.74 • 10 -4 


-4.23 ■ 10~ 5 


2.99 • 10 -4 


0.0149 


-5.90 


-5.32 


K2 max 


1.21 ■ 10~ 4 


3.77- 10~ 5 


-1.65 • 10~ 4 


-1.04- 10~ 4 


2.27- 10~ 4 


0.139 


6.18 


5.26 


K2 min 


9.03 ■ 10~ 4 


-1.19 ■ 10~ 3 


-2.15 • 10 -5 


4.56 • 10~ 4 


-4.49 ■ 10~ 4 


-0.556 


0.0338 


3.16 


K3 max 


-7.67- 10~ 4 


-1.23 ■ 10~ 3 


-8.03- 10~ 5 


-1.07- 10~ 4 


1.31 ■ 10~ 4 


-0.0348 


5.00 


6.42 


K3 min 


9.00 ■ 10~ 4 


-1.60 ■ 10~ 3 


2.30 ■ 10 -6 


-3.75- 10~ 4 


-5.64- 10" 4 


0.0656 


-0.0242 


0.0140 



TABLE IV: The same as Table 11111 for Model II. Now we also have the values of Ki as well as the extremal values of |k$|. 
The Ki are given in units of MeV. 
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tan 2 6*i2 


tan 2 623 


sin 2 6*i3 


Amli 


Am.23 


f -i m 'a y 

>L X IIIOjjY 


0.48 


0.80 


0.0061 


0.077 


2.4 




0.45 


0.79 


0.0012 


0.081 


2.5 


■f-2 llldA 


0.41 


1.19 


0.0007 


0.078 


2.5 




0.44 


0.93 


0.0031 


0.080 


2.4 




0.41 


1.06 


0003 


0.079 


2.4 


£i min 


0.47 


1.04 


0.0072 


0.081 


2.6 


■^2 IU33C 


0.44 


1.03 


OOOQ 


0.080 


2.7 


mm 


0.47 


1.04 


0.0075 


0.080 


2.6 


m&x 


0.44 


0.78 


0053 


0.079 


2.7 




0.44 


0.93 


0.0031 


0.080 


2.4 


k 1 in fi y 


0.40 


0.83 


0.0034 


0.080 


2.6 


M"i min 
r\j x iiiiii 


0.47 


0.81 


0056 


0.081 


2.7 


AC2 max 


0.49 


1.18 


0.0001 


0.078 


2.6 


K2 min 


0.42 


0.95 


0.0051 


0.081 


2.4 


max 


0.44 


0.74 


0.0014 


0.081 


2.8 


K3 min 


0.45 


0.69 


0.0055 


0.080 


2.6 



TABLE V: Model II (continued). Again, the values for 
Amli and Am 2 3 are given units 10 _3 eV 2 . 



[£' 3 max]. Thus we get the strongest constraints from 
the allowed mass ranges and from tan 2 #23- 

In the case of Model I, we only have five free pa- 
rameters. With these we must fit the two neutrino 
masses, two mixing angles and one upper bound. It 
is thus perhaps not surprising, that except for £' 1i the 
allowed ranges for the five parameters are quite nar- 



1 min 



is consistent with zero. In summary, we 



find from Table Un 



7.10 


10 


— 1 


< 




< 9.38 


10 


1.58 


10 


-3 


< 


\h\ 


< 1.78 


10 


1.00 


10 


-7 


< 


141 


< 9.27 


10 


3.61 


10 


-4 


< 


141 


< 4.28 


10 


5.04 


10 


-4 


< 


141 


< 5.71 


10 



(70) 
(71) 
(72) 
(73) 
(74) 



We shall employ the central values of these regions in 
the discussion of the resulting collider signals, below. 



2. Model II 



B. Discussion of the Results 

1. Model I 

In Table [nil we present the fit values for the param- 
eters £i, i'j in Model I. In the five columns on the right, 
we also include the resulting neutrino mass and mix- 
ing parameters. Of the large number of solutions we 
find, we present those where the parameters \£i\ take 
on extremal values. For example in the first row of Ta- 
ble HTU \£i I takes on the largest value we have found. 
Wc can now see in the five columns on the right, that 
tan 2 #12 and Amfi, are at the upper limit of their al- 
lowed ranges, Eqs. (JTJ) and (|64|) . respectively. This is 
as we would expect from Eq. (|6ip . where we see that 
£\ influences the first two generations. On the other 
hand, for example in the seventh row, where \£' 2 \ is 
maximal, we see that tan 2 #23 and Am§ 2 are at the 
upper limit of their allowed ranges. Similarly, in the 
second row, where \£\\ is minimal, we see that Am|i 
is at the lower end of its allowed range. For the eighth 
row, where \£' 2 \ is minimal, Am 32 is at the lower end of 
its allowed range. Overall, we see that tan 2 6*12 is at its 
upper limit for \£\ max] and also for [£2 min] . tan 2 #23 
is at its upper limit for [£' 2 max] and at its lower limit 
for [£' 2 min] and [£' 3 max] . sin 2 #13 is always well within 
its limits and thus does not pose a real constraint on 
our fit. However, we do predict a value between 0.003 
and the current upper bound. For Am 2 ^, we are at 
the upper end of the allowed range for [£i max], [£ 2 
max], and [£ 3 min]. We are close to the lower range 
for [£i min], [£[ max], and \£' x min]. For Ami, we are 
at the upper end of the allowed range for [£' 2 max] and 



In Table ITVl we show the results of our fit to Model 
II. We now have a total of eight free-parameters. We 
vary the values of the mass mixing parameters in the 
interval 

0.01 MeV < Ki < 10 MeV . (75) 

Due to the enhanced freedom, we see that we now 
have solutions, where are proportionality constants 
£i, i'j = 0(1O~ 7 ), which is consistent with zero, in 
our approach. We see that we push the upper bound- 
ary of tan 2 9i2 for \£\ max] and [K2 max] and the lower 
boundary for [ki max] . For tan 2 #23 we push the up- 
per and lower limits for [K2 max] and for [K3 min] , re- 
spectively From sin 2 $13 , we again have basically no 
constraint on our model, beyond those of the other pa- 
rameters, i.e. we are always well within the CHOOZ 
bound. 

We are pushing the upper end of the allowed range 
of Aniji for [£\ min,] \£\ min], [fti,2 min], and [K3 
max]. For Am^, we are at the upper end for [£' 2 
max], [£' 3 max], and [ki min]. We are at the lower 
end of Ato 2 ^ for \£\ max] , [£ 2 min] , [£[ max] , [£' 3 min] , 
and [K2 min]. Thus the mass ranges set the strictest 
limits on our parameters, the angles arc fairly easy to 
accommodate. 

Due to the enhanced freedom, we see that any one 
of our parameters can consistently be set to zero. This 
is of particular interest when trying to extract typical 
collider signatures, below. In Model II, it is thus diffi- 
cult to discern an identifying experimental signature. 
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Al22^211 


A132A3H 




A231A311 


Bound 


4.0 ■ 10~ 8 


4.0 ■ 10 -8 


4.0 ■ 10~ 8 


4.0 ■ 10 -8 





K12K21 


^113-^131 


X' X' 

A il3 A i31 


A 'lkl A 2kl 


A 'llj A 21j 


Bound 


10~ 9 


3- 10~ 8 


8- 10~ 8 


8.0- 10~ 8 


8.5 ■ 10~ 8 



TABLE VI: Bounds on the products of B3 couplings [611 ] . 
The first four and the last two bounds arise from contribu- 
tions to the process fj,Ti-^ eTi. The fifth bound arises from 
contributions to Attik and the sixth and seventh from con- 
tributions to Amj. 



VII. BOUNDS ON THE PRODUCTS OF THE 
PARAMETERS OF THE SIMPLE B3-MODEL 

Before discussing the consequences of our model, we 
first consider the low-energy constraints. Typical 
bounds on single B 3 couplings are of order 0.1 to 0.01 
[H, |70| . However, we necessarily have multiple cou- 
plings in our models and thus must take into account 
the bounds on products of couplings [5^, [6(| HH, , 
which are also typically much stricter, due to lepton 
flavor violating effects. The strictest product bounds 
(< 10" 7 ) of Table II in Ref. [H[ are given in TablcED 
We can now investigate whether our couplings sat- 
isfy the bounds in Table IVT1 For this we use the values 
in Tables HI and HT1 and assume tan/3 = 10. We find 



Al22A 211 


= 2.1 ■ lO" 6 ^ 


(76) 


^132 A 311 


= 


(77) 


Am A' m 


= -9.9 ■ 10 -9 € 2 ^'1 


(78) 


A 2 3lA 311 


= 


(79) 


X' X' 

A il2 A i21 


= 1.1-10~ 7 «) 2 


(80) 


X' X' 
A 113 A 131 


= 2.7- l0- 9 (£[) 2 


(81) 


X' X' 


= 2.7- 10" 9 (4) 2 


(82) 


X' X' 


= 1.2 • lO^ 7 ^, k = 1 


(83) 


A W A 21j 


= 2.0- 10" 6 ^ 2 , j = 2 


(84) 



In the last two cases the chosen indices result in the 
largest possible value. The strictest bounds result 
from Eqs. (HD and $EU§: i x l' 2 < 0.02, and l\ < 0.1. 
From Tables IITlI and ITVl we see that these are always 
satisfied in our numerical solutions. 



VIII. COLLIDER TESTS 

An essential feature of B3 neutrino models, is that 
they necessarily lead to observable consequences at 
colliders. Resonant slepton or squark production re- 
quires couplings of order 10~ 3 or larger [57] . As can be 
seen from Tables HI and [TTI together with the numerical 
results presented in Tables IIII1 IIV1 this is not possible 
in our models. However, it is well known, that at the 
LHC squark and gluino production provide the largest 




FIG. 2: The two-body decay of a scalar particle into two 
fcrmions. 



supcrsymmetric cross sections. This is independent of 
whether or B3 is the relevant symmetry. The pro- 
duced squarks and gluinos then cascade decay within 
the detector to the LSP. In particular, this also holds 
for the BC benchmark points (58[, where fi is the 
LSP. In this paper, we focus on stau-LSP scenarios, 
as outlined in Ref. [58[ . We shall focus on the essential 
features, a full phenomenological analysis goes beyond 
the scope of this paper and will be presented elsewhere 
[7ll ]. The neutralino LSP case requires a full treat- 
ment of the scalar potential in order to determine the 
relevant couplings and masses. This in turn requires 
assumptions about the soft-supersymmetry breaking 
sector, which also goes well beyond the scope of this 
paper. We shall consider this elsewhere (7l| . 



A. Stau LSP Decays 

As discussed in detail in Refs. [Tll.[5c|. there are ex- 
tensive regions of mSUGRA parameter space, where 
the scalar tau is the LSP. The final state collider sig- 
nals will be determined by the dominant decays of the 
stau. The lightest stau, fi, is an admixture of right 
and left stau. 



ri 



cos t)f tr + sin 6f Ti 



(85) 



with Of the mixing angle. In Ref. [58j, it was found 
that in the representative benchmark points (BC1- 
BC4) the fi is dominantly a right-handed stau with 
\9f I < 0.3 (in radiants), i.e. the fi-LSP is more than 
91% t r and sin 2 6? < 0.09. 

In our model, we have a wide range of non-zero B3 
couplings, where the corresponding operators couple 
directly to the stau. The stau can thus decay via the 
two-body mode into two spin-1/2 fcrmions /i ; 2 shown 
in Fig. O The correspon ding partial decay rate is in 
given in leading order by [72[ 



r(n^/x/ 2 ) - 



jv c |A| 2 e 2 p c 

8ttM? 

A^ c |A| 2 e 2 

16tt 



■{Ml -m\ -ml) (86) 



(87) 



where 



P 2 

rem 



[Ml - (mi + m 2 ) 2 ] [Ml - (m x - m 2 ) 2 ] 



4M? 
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Operator 


Decay Mode 


rffr -> bc)/M f , 


L\ L$ E\ 


— > e~ u e 


1.5 • 10" 12 ^| 


L\ L'$ E'$ 


tT — > T — v e 


2.0 • 10~ 4 


L\ L% E$ 


— > e~ v T 


1.8 • 10" 4 


L2 L% E2 


t7 —> uTv,, 

1 r* l-i 


6.7 • 10~ s £\ 


L2-L3E3 


T7~ — > T~V a 

1 H- 


2.0 ■ 10~ 4 £% 


L2L3E3 


ff — > /i~^ T 


1.8 • 10~ 4 £| 


L S Q 2 D 3 


fj~ — > cb 


5.4- 10" 7 (^) 2 


L3Q2D2 


f {" — > c s 


1.8- lO" 7 (^) 2 


UQ X D 2 


-ms 


1.1 • 10" 8 {£' 3 ) 2 



TABLE VII: Decay modes and partial decay widths of a 
stau LSP given as a function of the relevant li, i'j. In the 
second and the fifth decay modes we have added the two 
contributions from the doublet and the singlet stau. We 
have set 8f = 0.3 (in radiants) as obtained for BC1. 



mi,2 denote the final state masses. A denotes one of 
the following B3 couplings relevant for tree-level stau 
decay 



A G {A131, A133, A232, A233, A 3 fe } , 



(89) 



and N c is the colour factor. N c = 1 for the decay via 
the LLE operators, and N c = 3 for the LQD opera- 
tors. G = (cos 6f, sin Of), depending on whether the 
?i couples via the right- or the left-handed stau com- 
ponent. Eq. (|87[) is taken for the case where mi. 2 — > 0. 
This is a good approximation for Mf x < m t0 p, which 
is the case for all the BC benchmark points. 

Given the above decay formula we can now compute 
the decay rates for the dominant decay modes using 
the numerical values in Tables [T] and |TT] for the rele- 
vant coupling. We expect the decays where the right- 
handed stau component couples directly to dominate, 
due to the small mixing angle in the stau sector. Fur- 
thermore, for A333 which is a potentially large coupling 
the large top quark mass kinematically blocks the de- 
cay, for the stau masses we consider here. We present 
the results for the decays in terms of the li, l'^ in Table 
IVIII For completeness, we have included the couplings 
involving £3 , which we have neglected in our neutrino 
parameter fits. We see that for substantial decays via 
the corresponding operators, we would require, e.g. 

h » k,2- 



1. Model I 



If we consider the BC1 benchmark point, we have 
M fl = 148.38 GcV. Using Tables HH [TV] we can then 
compute explicit values for the partial widths and the 
branching ratios. In Model I, |4| « 8 ■ 1(T 4 , \£ 2 \ « 
2-10 -3 and \£' 3 \ « 5-1CP 4 . Using these values, we then 
obtain the total decay width and lifetime in Model I 
at BC1 

r(fi) = 260 eV, T(n) = 2.5 • 10~ 18 scc. (90) 



We see that the stau-LSP always decays within the 
detector. It also will not lead to a detached vertex. 
For the branching ratios of the decay modes in Table 
IVIII we obtain 



Br(fi - 


-> r" 


~v e ) 


= 0.072 


(91) 


Br(fi - 


-> 


~v T ) 


= 0.065 


(92) 


Br(fi - 


-> 




= 0.45 


(93) 


Br(fi - 




~V T ) 


= 0.41 


(94) 



The other decay modes are negligible. The branch- 
ing ratios are independent of Mf t , i. e. in Model I 
they only depend on the £i, £'j. As the neutrinos are 
not visible, we have a combined branching ratio into 
charged tau leptons of about 52%. For squark or 
gluino pair production, we expect two stau's in the 
decay chains. The probability for then having two 
charged electrons/muons in the final state is about 
23%. Since the gluino/squark pair production cross 
section is very large for accessible supersymmetric 
masses, this should lead to an easily visible signal rate. 



2. Model II 

In Model II, we have eight free parameters and thus 
a much larger freedom. Furthermore for K3 7^ the 
stau can mix with the charged Higgs boson, leading 
to additional decay modes. In order to compute these 
properly, we must minimize the full scalar potential. 
This is beyond the scope of this paper. We can es- 
timate the stau-Higgs mixing to be k^/^jl. Using the 
Feynman rules in Fig. 8 of Ref. [73[ for the charged 
Higgs coupling to the tau lepton, we then typically 
find a product of mixing times couplings of order 10~ 7 , 
for « 3 = 1 MeV, tan/3 = 10 and (i = 200 GeV. This 
would lead to an additional decay f — ► rv. The cou- 
plings to the second generation quarks are another 
order of magnitude smaller. These Higgs-mixing cou- 
plings are negligible compared to the direct stau de- 
cay couplings, in most cases. However, in general they 
must be included. A proper complete treatment will 
be given in Ref. [tT| - Here it shall suffice to present 
one example case from Table IIVI employing only the 
direct decays from Table IVHl We choose the example, 
such that <C £' 3 , which differs from Model I. 

We consider the case: [£[ min], where £\ = —4.65 • 
10~ 6 , £ 2 = 9.62 • 10~ 6 and £' 3 = -4.04 • 10~ 4 . We find 
for the total decay rate and the lifetime 



r(n) = 2.4-10" 
r(fi) = 2.7-10" 



; eV. 



sec . 



(95) 
(96) 



We see that the width is now substantially smaller 
and thus the lifetime correspondingly larger. For a 
Lorentz boost 7^ = 10, we have a decay length of 
about 100/im. This is on the boarderline of visibility 
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for a detached vertex, 
find 



For the branching ratios we Higgs Yukawa couplings via a small set of parameters 

£i, £' . This results in simple relations between the B3 



Br(fi -»■ hadrons) = 0.73 (97) 

Br(fi -> r~u e ) = 0.026 (98) 

Br(fi -> e~v T ) = 0.023 (99) 

Br(fi -> r - ^) = 0.11 (100) 

Br(fi -> /ti _ i/ T ) = 0.10 (101) 

The stau now dominantly decays hadronically. In 
this specific case, we can have still a roughly 25% 
branching ratio to charged leptons. Or a probabil- 
ity of roughly 6% for two charged leptons in the final 
state. However, recall that we only scanned couplings 
down to 10" 7 . Thus we would expect solutions with 
even smaller £ 1<2 . In this case wc would have purely 
hadronic final states and we must resort to the tech- 
niques used in Ref. [74|, where the UDD R-parity 
violating operators were studied. 



IX. CONCLUSIONS AND OUTLOOK 

We have presented a simple ansatz for the B3 Yukawa 
couplings, relating them directly to the corresponding 



couplings presented in Tables HI and HT1 We have given 
estimates of these parameters in order to obtain the 
correct neutrino masses and have then numerically de- 
termined the precise values. These are summarised in 
Tables IIH| ITVl Wc then discussed the resulting col- 
lider signals for the case of a stau LSP. Depending on 
the fit values, we have found a wide range for the pos- 
sible branching ratios of the stau-LSP. In forthcoming 
work, we shall give a detailed investigation of how to 
disentangle these models at the LHC. 
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